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Abstract 

CoUisionless regime kinetic models for coherent nonlinear Alfven wave dynam- 
ics are studied using fluid moment equations with an approximate closure anzatz. 
Resonant particle effects are modelled by incorporating an additional term repre- 
senting dissipation akin to parallel heat conduction. Unlike coUisional dissipation, 
parallel heat conduction is presented by an integral operator. The modified deriva- 
tive nonlinear Schrodinger equation thus has a spatially nonlocal nonlinear term 
describing the long-time evolution of the envelope of parallel-propagating Alfven 
waves, as well. Coefficients in the nonlinear terms are free of the (1 — (3)~^ singu- 
larity usually encountered in previous analyses, and have very a simple form which 
clarifies the physical processes governing the large amplitude Alfvenic nonlinear dy- 
namics. The nonlinearity appears via coupling of an Alfvenic mode to a kinetic 
ion-acoustic mode. Damping of the nonlinear Alfven wave appears via strong Lan- 
dau damping of the ion-acoustic wave when the electron-to-ion temperature ratio 
is close to unity. For a (slightly) obliquely propagating wave, there are finite Lar- 
mor radius corrections in the dynamical equation. This effect depends on the angle 
of wave propagation relative to Bq and vanishes for the limit of strictly parallel 
propagation. Explicit magnetic perturbation envelope equations amenable to fur- 
ther analysis and numerical solution are obtained. Implications of these models for 
coUisionless shock dynamics are discussed. 
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I. Introduction 

It is widely believed that Alfven waves play an important role in in- 
terplanetary plasmas. High level of magnetohydrodynamic (MHD) wave 
activity and shock waves are observed in the planetary and solar wind 
plasmas [1-3]. The plasma environment is characterized by a quite weak 
mean interplanetary magnetic field and large amplitude magnetic field 
fiuctuations caused by Alfven and magnetosonic waves. Thus, nonlinear 
effects play an important role in the evolution of such waves. Interest in 
large-amplitude Alfven wave dynamics arose from the attempts to under- 
stand the steepening of wave trains and shock formation in space as well 
as from general interest in nonlinear waves as a whole. 

Previous studies have shown [4-6] that nonlinear Alfven and fast mag- 
netosonic waves are described by the derivative nonlinear Schrodinger 
equation (DNLS). In this model, a parallel ponderomotive force of a high- 
amplitude wave perturbs the plasma density. The force squeezes plasma 
out the regions of larger magnetic field, locally decreasing the plasma 
density, and increasing the Alfven wave velocity. As the amplitude of a 
circularly polarized Alfven wave varies on time scales much slower than 
the wave frequency, the section of a wave train with higher amplitude 
propagates faster than the part with lower one. This gives rise to nonlin- 
ear steepening of the large-amplitude Alfven wave. Dispersion, as usual, 
ultimately controls steepening. The DNLS thus describes the long time- 
scale dynamics of the envelope of such Alfven waves. The DNLS is an 
integrable equation, which describes soliton and multi-soliton solutions, 
shock waves, modulation instability of solitons [7-9] etc. As was shown 
by Longtin and Sonnerup [10] and Wong and Goldstein [11] (from anal- 
yses of the full set of MHD equations), modulation instability occurvs 
for left-hand polarized waves if /3 < 1 and for right-hand polarized waves 
if P > 1. This DNLS modulation instability is strongly sensitive to the 
sign of the coefficient of the cubic nonlinearity, which is propotional to 
(1— Indeed, at /3 = 1 (typical of solar wind plasma), the sound speed 
approaches the Alfven velocity, and resonant energy exchange between 
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sound and Alfven waves is strong, in the absence of parallel dissipation. 
Then, the Alfven wave is always in phase with the acoustic wave, giving 
rise to fast steepening of the front of the nonlinear Alfven wave train. The 
rate of the steepening is roughly proportional to {l—/3)~^. Significant con- 
version of Alfven wave energy into ion thermal energy (i.e. ion heating) 
occurs as well. It is interesting to note that the overwhelming prepon- 
derance of theoretical work in this field of quasi-parallel Alfvenic shocks 
and solitons is based upon a simple MHD plasma model. Thus, with 
the notable exceptions of Ref. [6,12-14], the theoretical 'lore' of nonlinear 
Alfven waves in a coUisionless, /3 ~ 1 plasma is built upon a conceptual 
paradigm constructed for a collisional, (3 < 1 system. This state of affairs 
is due, in part, to the intractability and unwieldiness of straightforward 
kinetic analysis of nonlinear Alfven waves. In this paper, we offer a sim- 
plier approach which exploits recent developments in the theory of fluid 
modelling of kinetic effects. 

We can expect that the influence of kinetic effects, such as Landau 
damping must become important when /3 approaches unity. A fully ki- 
netic calculation was undertaken by Register [5] for the case of a high-beta 
plasma. For parallel propagation, his results counsidc with equations ob- 
tained later by Mj0lhus and Wyller [12, 6] and Spangler [13, 14]. The 
effect of Landau damping appears in the DNLS via an additional cubic 
term which is an integral operator over space. There is no (1 — singu- 
larity in the nonlinear derivative term. The nonlocality is a consequence 
of the finite time history of an ion transit through an envelope modulation 
(i.e. the coherence time of an ion with the modulation is not infinitesimal 
in comparison to the modulation growth time). Thus nonlinear Alfven 
wave dynamics are nonlocal, and not governed by local spatial derivatives 
of the perturbed field, alone. The other attempts to include Landau 
damping in the nonlinear Alfven wave evolution used kinetics to calcu- 
late second-order density (or pressure) perturbation, but the DNLS-like 
equation (for b) was obtained from MHD equations. Mj0lhus and Wyller 
[12, 6] used a guiding center formalism and Spangler [13, 14] used the 
full Vlasov equation. Both methods predict a nonlocal integral term, and 
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the coefficients of the nonhnear terms coinside with those of Rogister [5] . 
In these cases, the coefficient of the derivative cubic nonhnear term does 
not change sign at = 1, except for large values of T^/Ti (the electron- 
to-ion temperature ratio). The coefficient of the nonlocal term is always 
negative. Both coefficients of the nonlinear terms (cubic and nonlocal) 
depend strongly on j3 and on T^/Ti. 

Obviously, a simple fluid model of collisionless shocks is desirable for 
reasons of tractibility. However, kinetic effects are essential for describing 
interplanetary plasma dynamics. It is known that fluid models poorly 
represent most kinetic effects. Several authors have suggested that kinetic 
effects, such as Landau damping may be modelled in fluid equations by 
adding parallel dissipation terms [15,16]. Recently, a closure method for 
modelling of kinetic effects was developed [17-21] and applied to ion- 
temperature gradient instabilities using gyrofiuid models. This closure 
anzatz for fluid moment equations i) ensures particle, momentum and 
energy conservation, ii) takes a simple form in wave-number space, and 
iii) has a linear response function very close to that of a collisionless, 
Maxwellian plasma. In this paper we use the simplest version of the 
method, namely that proposed by Hammet and Perkins [17]. 

In this paper we investigate effects of dissipation on the nonlinear, 
parallel-propagating Alfven wave. We first use one-fluid MHD equations 
with parallel dissipation coefficients //y and x\\ which are constant and 
independent of mode frequency and wave-number. Later we replace con- 
stant x\\ with the integral operator representation of Hammet and Perkins 
[17], and obtain a modifield DNLS similar to that obtained by Mj0lhus 
and Wyller [12, 6] and Spangler [13, 14]. However, our approach yields 
expressions for coefficients of the nonlinear terms which are much simpler 
than theirs and allows clear, unambiguous physical interpretation of the 
results. We also derive the Te/T^ dependence using two-fluid MHD equa- 
tions. In an ion-acoustic wave, ion density perturbations of the electron 
background are mediated by electric field effects. In collisionless plasmas 
the ion-acoustic branch, thus, replaces the acoustic branch as caused by 
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gas-kinetic pressure perturbations. The results indicate that the nonlin- 
earity of Alfven waves in colHsionless plasma is controlled by the coupling 
of the Alfven mode to an ion-acoustic mode, i.e. an ion-density pertur- 
bation mediated by an electron responce electric field. This is in contrast 
to the conventional view [4,7,8] of this nonlinearity as due to a pondero- 
motive (i.e. a gas-kinetic/collisional) plasma density perturbation. The 
kinetic damping on resonant particles is nothing more than the usual 
strong Landau damping of an ion-acoustic wave for Te — Tj. This damp- 
ing leads to enhanced ion heating which further raises Tj/Tg. 

We also consider the case of a slightly obliquely propagating wave. 
There, one can expect that other kinetic effects may be relevant to the 
dynamics. It is shown that the E x B drift in the electric field produced by 
charge separation in an ion-acoustic wave with an ambient magnetic field 
is not significant. Another effect which is important, however, is caused 
by gyro-averaged electric and magnetic fields acting on a particle over 
the scale of a Larmor orbit. (There are no {k±pi)'^ effects on dispersion 
in the DNLS approximation. Such effects may appear in higher-order in 
B/Bq calculations.) This effect enters the modified DNLS in a way similar 
to coUisional dissipation. The finite Larmor radius correction depends on 
the angle of propagation of a wave and disappears for a strictly parallel- 
propagating waves. We explore modulational stability for the general 
case of a dissipative nonlinear Alfven wave. This instability is important 
when one considers the origin and evolution of solitons, wave packets, 
shock waves, etc. 

The rest of this paper is organized as follows. In Section II we de- 
rive the evolution equation for dissipative nonlinear parallel-propagating 
Alfven waves. In Section III the modified DNLS with the resonant par- 
ticle effect is obtained. In Section IV we consider the influence of finite 
Larmor radius corrections to the dynamics of a slightly oblique, nonlin- 
ear Alfven wave. In Section V we investigate the modulation instability 
of dissipative nonlinear Alfven waves. Section VI is a discussion of the 
results obtained. 
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II. Dissipative Nonlinear Alfven Waves 



We obtain the equation governing nonlinear wave dynamics from a 
multiple time scale expansion of the dissipative MHD equations. Dissi- 
pation is included via the parallel viscosity //|| or the parallel heat con- 
ductivity x\\- We considered both these cases for the following reason. 
The first case, /i 0, corresponds to the three-moment fluid model of 
Landau damping. The four-moment fluid model should contain both x\\ 
and to close the equations. As shown in Ref. [17], the best fit of 
the linear-response function of of this model to the linear-response func- 
tion of a Maxwellian plasma is achieved when //y = and X|| 0- This 
is the second case we consider here. Note that any closure of this sort 
which tacitly assumes a Maxwellian plasma intrinsically fails to capture 
the strong local modification of the distribution function associated with 
large-amplitude turbulence, trapping, etc. 

In the derivation, we follow Ref. [4]. We assume a plane wave prop- 
agating in the z direction. Thus, all quantities are functions of z and t. 
The MHD equations are written as: 

dp du dp ^ 

du du dp dS d h'^ ^ 

dv dv Bn dh ^ 
dh dh du „ dv 

where the velocity components are v = V±,u = V\\ , the magnetic field 
components are b = Bx, Bq = B\\, p is the mass density, p is the pressure, 7 
is the polytropic constant. Also, 5" and q are the parallel momentum and 
heat fluxes respectively, i.e. S = —pp\\{du/dz) and q = —nx\\idT/dz), where 
n = p/rrii and T is the temperature. There is no longitudinal perturba- 
tion of the magnetic field, as V • B = 0. Eqs. (l)-(5) are the equations of 
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continuity, longitudinal momentum, energy conservation (with q = be- 
ing the adiabatic equation), transverse momentum and transverse flux, 
respectively. Dissipation, as stated above, appears only via either mo- 
mentum flux ^ 0,x = 0) or heat flux (fj, = 0,x We expand the 
system of Eqs. (l)-(5) in powers of e = b/Bo. To avoid secularities at third 
order, we assume multiple time-scale dependence, i.e. t(2n) ^ te'^^ are the 
independent variables. 

The zeroth-order solution of the system (l)-(5) is the equilibrium: 
Po = const , Vo = 0, Bq = Boeg , where is the unit vector in z direction. 
The flrst-order solution of Eqs. (4) and (5) is the linearized Alfven wave, 
i.e. bl = bi(2;±v^t(o),t(2)), vi = (vA/-Bo)bi + vi(t(2)) , where v\ = BI/A-kpq is 
the Alfven velocity. Eqs. (l)-(3) also have dissipative acoustic waves as a 
solution. As we consider travelling hydromagnetic waves propagating to 
the right, we set pi = ui = and and choose a minus sign in the argument 
of bl. 

At second order, we consider py ^ and X|| 7^ models, separately, 
a) p|| ^ model. 

Eqs. (l)-(3) become, in second order: 

dp2 ^ du2 , . 

du2 , dp2 d'^U2 d hi 

{P2 - cIp2) = 0, (8) 



^^(0) 

where = 7Po/Po is the sound speed. Prom this system, excluding the 
free sound wave solution (as in flrst order), we have: 

U2=—p2 + U2{tt2)), (9) 

Po 

P2 = ^^ip2 + ^||^'A^ - (10) 
P2 = C^P2+P2(i(2))- (11) 
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dz', (12) 



Solving this system of equations, we, finally, arrive at the equation for 
the second-order velocity perturbation: 

where /3 = (?a/v\ and = [^\\/va){^ — ■ This is the inhomogeneous 
solution of Eq. (11), finite for all —00 < z < +00. The term (6f ) appears in 
Eq. (12) from density conservation: d{p) /dt = 0. Note that it is the devi- 
ation from the mean ponderomotive force which causes density bunching 
and thus wave steepening. Thus (p) = po, so that (142) = {va/ Po) (P2) = 0. 

b) X|| 7^ model. 

Eqs. (l)-(3) in this case are replaced by: 

dp2 du2 , . 

;^^+Po-^ = 0, (6) 
at(o) dz 

^°at(o) dz dzSir' 

{p2-CaP2) = - (8') 



dt^o) dz' 
The conductive heat flux can be represented as follows: 

92 = -noxii - mxii ^Ti - n2X|| ^Tq. (13) 

The second term here vanishes because ni = pi/rrii = 0, as we set in the 
first-order equations. The third term vanishes also because we consider 
only a homogeneous plasma temperature Tq = const. The second-order 
temperature pertubation is, in turn, given by: 

no no 

where = To/rrii is the thermal velocity of the particles. Now, from 
Eqs. (6), (7'), (8'), and using (13) and (14) we obtain: 

W2 = — P2 + ^2(^(2)), (9) 
Po 

P2=v\p2 - ^+P2{t(2)). (10') 

{P2 - cIp2) + ^|-(P2 - ^P2) = 0. (11') 

Va oz 7 
8 



Prom these equations we obtain the second-order velocity perturbation 
for the X|| 7^ model as: 

""'^^^ 2(1-/3/7)1 B', 



+ ^(1 



XII 



Here = (x||/t;A)(l — /3/7)/(l — /3)- The velocities W2(^) and 'U2(x) are 
expressed in the frame co-moving with a wave (i.e. z z — VAt(Q))- The 
coefficient 7 = 3 is chosen to ensure energy conservation [17]. 

Eqs. (4) and (5) are, at second order, equations for linear Alfven 
waves. One may easily show that for all t(2) one may set b2 and V2 to zero 
[4]. At third order, upon substituting first- and second-order solutions 
into Eqs. (4) and (5) one obtains the equation for dhi/dt(^2)- 

1 a 

For a dispersive term [6,24] one must invoke finite Larmor radius effects. 
Since Eq. (15) describes the transverse components of b, there is no in- 
fluence of parallel dissipation effects on the linear dispersive term. We 
thus add a general form of this term to obtain a DNLS-like magnetic per- 
turbation envelope equation. Introducing r = t(2) and (f) = {bi^ + ibiy)/BQ 
we have: 

where Qi = eBo/rriiC is the ion-cyclotron frequency, upper (— ) and lower 
(+) signs on the third term refer to right and left elliptically (circularly) 
polarized waves. Eq. (16) with (12) or (12') inserted for ui describes the 
coherent nonlinear dynamics of an Alfven wave train subject to damping. 
Damping results in the appearance of a new nonlocal integral contribution 
to the envelope equation. We define the integral dissipation operator as: 

Jl[F]{x) ^e-""/^ J e'''/^F{x')dx', (17) 
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where F is the arbitrary function the operator acts on and L is a char- 
acteristic length proportional to the dissipation coefficient (i.e. L ~ 
or L~x||/b^(l-/?)] )• 

From comparison of Eqs. (12), (12') to the result of Ref. [13] we con- 
clude that three-moment approximation {x\\ 7^ model) of Landau damp- 
ing recovers the correct functional dependence U2 — c\h\ + C2>C[&i], where 
£ is a nonlocal operator, and ci and C2 are some coefficients. However, 
the integral operator in this case (Eq. (17)) is different from the resonant 
particle operator of Ref. [13] (see Eq. (21) in the next section). The two- 
moment approximation (//|| 7^ model) describes the nonlinear physics 
incorrectly. Indeed, it lacks a free cubic nonlinear contribution to the 
expression for U2. This is in agreement with the conclusions of Ref. [17], 
which, stated simply, are the higher the moment approximation used, 
the better the description of Landau damping in a fluid model which 
results. In particular, the linear response function of the three-moment 
fluid model is much closer to the exact Maxwellian linear response func- 
tion then that of the two-moment fluid model. From now on we use the 
three-moment fluid model (xy ^ 0). 



III. Coherent Nonlinear Alfven Waves with Landau Damping 

Using the X 7^ model, kinetic effects are modeled by a longitudinal 
heat flux. In the linear closure approximation, this flux is (in wave- 
number space [17]): 

92fc = -nQXi^^^^ikT2^. (18) 

Here the temperature perturbation is given by Eq. (14), and xi is a di- 
mentionless fit coefficient for the model. The choice xi = 2/\/7r gives the 
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best fit to linear Landau damping. By performing the inverse Fourier 
transform of q2i^, we obtain the real-space representation of q2 as: 

q2(z) = ^= lim / dke-\^\^^^^q2, 

noXiV2vt rT2{z + z')-T2iz-z') 



TT 



dz' 



-^T2{z')dz'. (19) 



./-oo Z' 



Here V stands for the Cauchy principal value integral. We have also added 
the factor of exp(— |A;|5) to control the otherwise infinite integral. Prom 
comparison of this equation with Eq. (13), we conclude that the effect of 
resonant particles in our theory is reduced to the following replacement: 

^|--X||>C. (20) 
Va oz " 

We thus define here the resonant particle (nonlocal) integral operator: 

C[F]{x)=- F{x')dx' (21) 



^ •'— oo ^ 



with the coefficient = Xiv^vt/^A • 

We now consider properties of the resonant particle operator in more 
detail. We need to know the inverse operator, so that = 

1. It is easy to find in /c-space representation, i.e. C~^C = C^^ik/\k\ = 
C^^i sign(/c) = 1, thus, C^^ = —i sign(/c) and, returning to real space: 

= -C. (22) 

Formula (22) reflects the fact that Landau damping is time reversable, 
i.e. that a system returns to its initial state when time is reversed. There 
is no microscopic information loss in a Vlasov system, unlike in the case 
of collisional damping. Writing an evolution equation in the form {d/dt + 
C)'4^ = 0, we see that after time inversion {t —>■ —t), C must satisfy Eq. (22) 



11 



to represent damping. It can also be shown that the operator L commutes 
with any standard differential or integral operator, i.e. 



d_ 

dz 



d_ 



F{z) 



{C[F{z)])dz = C I F{z)dz 
C[c] = 0, 



(23) 



where c = const. 



Making a replacement, Eq. (20), in the energy conservation equa- 
tion (11'), we have: 



H - {hi) 



(24) 



This equation can be easily solved for p2 using Eq. (22). The second-order 
velocity perturbation is then: 



l$l 



]}' 



(25) 



U2 = va|Mi |$|^ + M2C 
where |$|^ = - = (fcf _ {hl))/Bl, and the factors Mi and M2 are: 

(26a) 
(266) 



,,(1) 1 (l-/^)+xg(l-/?/7) 

1 2 (1-/3)2 + £2(1 _^/^)2' 



7-1 



7 (l-/3)2 + xg (1-/3/7)2' 



Here xy = xiV^^^Jli 7 = 3 and xi = according to Ref. [17]. The 

superscript (1) means that the result refers to the one-fluid model. 

Eq. (16) together with Eq. (25) recovers the modified derivative non- 
linear Schrodinger equation governing the envelope of nonlinear Alfven 
modes subject to Landau damping. Fig. 1 represents the dependence of 
the factors M^^"^ and on 1//3 for comparison with the work by Span- 

gler [13]. This case corresponds to the T^/Ti = 1 case of Spangler. The 
coefficient Mi of Spangler is positive, always. The coefficient Mi^^ in our 
analysis changes sign at /3 = 1, although where negative, its absolute value 
is small. The disagreement is related to the lack of higher moments in 



12 



the closure scheme we use. The coefficients M^'^ and M2 of Spangler are 
negative and look alike. Analytical expressions for the coefficients in our 
approach are much simpler than those obtained from kinetic calculations. 
The phenomenon which is described here by Eqs. (25), (26) is typical for 
other resonance phenomena. The sharp resonance without damping (here 
the resonance between Alfven and sound waves) is, of course, smoothed 
by increasing damping. In other words, there is resonance broadening 
caused by the interaction of particles with a wave. In this respect, the 
coupling to dissipation calculated here is the coherent analogue of non- 
linear ion Landau damping. 

Until now we have considered only one-fluid magnetohydrodynam- 
ics. Based on this model, we derived the Alfven-sound wave resonance 
broadening in a DNLS-like equation. Astrohysical plasmas are essen- 
tially collisionless, i.e. the wave-length is much larger than mean free 
path, A <^ imfp- Thus, an ion-acoustic mode replaces the sound wave in 
the collisionless regime. Hence, it is natural to adopt the perspective that 
the nonlinearity of finite-amplitude Alfven waves arises from an ion den- 
sity perturbation, that is from coupling of the ion-acoustic and Alfven 
branches of the plasma oscillation. In other words, we must consider 
electrostatic density perturbations instead of pure neutral density per- 
turbation as in a sound wave. The electric fields of the ion-acoustic wave 
play a role in the dynamics of oblique nonlinear Alfven waves, as will be 
shown later. To include the ion-acoustic wave coupling effects instead 
of the gas-kinetic pressure perturbation, we must use two-fluid dynamics 
for both electron and ion fluids. The longitudinal momentum equation, 
Eq. (2), for ions in a collisionless plasma is: 

du du dip 9 ^ ,^„. 

where <^ is the electric potential. We neglect wave damping by electrons, 
which is always much smaller than damping on ions, i.e. Select. I lion — 
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mf./mi. Due to the quasineutrality of ion-acoustic waves, the first nonva- 
nishing term of the expansion of the Boltzman distribution for electrons 
is: 

ni ~ ne = no — . (28) 

e 

Upon substituting this into Eq. (27), and formally defining (for ions) 
P2 = n2Toe = (Te/Tj)p2, we can write Eq. (7') as: 

^ du^,Tedp2__d_bl 

Multiplying Eq. (8') by the factor T^/Ti and rewriting it with P2 we have 
for this new pressure the same equation, with the c^^ = {Tf,/Ti)c^g, - the 
speed of an ion-acoustic wave. We do the same with the closure equation 
(14). Thus, all equations still apply for the new pressure P2 when we 
substitute 

- (30) 

everywhere except in the X|| coefficient. There j3 appears as a combination 
of the ion thermal velocity and the Alfven velocity. 

When TpjTi ^ 1 an ion-acoustic wave experiences almost no damping 
(neglecting the feeble damping by electrons). In the opposite case, Tg/Tj ~ 
1, there is strong ion Landau damping. This is the case of the one-fluid 
calculation. So, we must redefine the fit coefficient: 

XI -XI /(I) (31) 

with some correction function /(Tg/T^) that reflects the behavior in two 
asymptotic regimes. This function must reflect the temperature depen- 
dence of the damping rate, and must equal unity at Tg/Tj — 1. The 
damping rate of ion-acoustic waves in a Maxwellian plasma is well known 
[22]: 

Then, for the correction function / we obtain: 



/( 
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Upon substitutions given by Eqs. (30) and (31) in Eqs. (26), we arrive 
at: 

(i-ft/s) +JS(i-ftf) 

= - Ix, S^f ^ 7 r^^— -2 (346) 

with X|| = Xi\/W7l f{Te/Ti). Figs. 2a and 2b show the functions Mi and 
M2, respectively, plotted vs. 1//3 for different electron-to-ion tempera- 
ture ratios {Te/Ti = 1,3,5,10). These figures closely resemble the graphs 
of Spangler [14] and Mj0lhus and Wyller [6], however the expressions 
Eqs. (34) are much simpler than their counterparts in Refs. [6,14], and 
easily amenable to detailed calculation and evaluation, etc.. The coeffi- 
cient M2 is always negative and has a peak near /?o = Ti/T^. There is slight 
difference when ~ Tj. The coefficient Mi should be positive in this case, 
for almost all values of (3. This difference can result from insufficient accu- 
racy of the three-moment approximation {x\\ 7^ 0) when Landau damping 
is very strong. To get a better fit one should take higher-moment approx- 
imations. Of course, for large amplitude Alfven waves the assumption of 
a Maxwellian plasma fails, as well. 

IV. Slightly Oblique Nonlinear Alfven Waves with Landau Damping 
and Finite Larmor Radius Effects 

In the previous sections we considered Alfven waves propagating 
strictly parallel to the magnetic field lines. Here we try to generalize 
our approach to waves propagating at a small angle to the magnetic 
field. It has been shown [12, 23-25] that dissipationless oblique nonlinear 
Alfven waves are governed by the DNLS equation, as well. We investi- 
gate here the influence of kinetic effects on wave dynamics. To include 
kinetic effects in the fluid model of the DNLS we use a two-fluid gyrofiuid 
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model [18,19]. Here we construct the simplest model which includes ki- 
netic effects such as finite gyro-frequency, Landau damping, etc. We 
omit gyro- viscous corrections to the Reynolds stress. These corrections 
are small for the slightly oblique propagation case in straight field ge- 
ometry, and can be incorporated (approximately) into the viscous model 
considered in Section II. 

Typically, the gyrofiuid equations are written in components pro- 
jected onto the magnetic field direction, 2;', and perpendicular to it. It 
is more convenient to us to work in the frame of a wave propagating in 
the z direction, which makes some angle <^ 1 relative to the z' direc- 
tion. We choose the y direction coincident with y' . We leave variables 
'unprimed' when measured in the frame of a wave, and use primes for 
variables measured in the frame of the ambient magnetic field. We con- 
sider contributions of kinetic effects separately for every equation of the 
system (l)-(5). The equation for the transverse flux is obtained from the 
Maxwell equations V x E = (-l/c)aB/at, V x B = (47r/c)J and the Ohm's 
law. We write the generalized Ohm's law as follows: 

1 777 ■ 777 ■ 

E = — [v X B] + ^[J X B] - ^Vp + V^ei- (35) 
c ecp zep 

There are no contributions from the last two terms. There is no contribu- 
tion from the term V</? x Vp because <^ is the first-order perturbation, and 
we have set po = const, pi = 0. The first term has already been accounted 
for. The second term leads to an additional term on the right-hand-side 
of Eq. (5). We write it in components as follows: 

The terms on the right-hand-side give rise to a linear dispersion term in 
the DNLS. This term describes dispersion due to finite ion Larmor radius. 
In gyrofiuid models there is no correction to the transverse momentum 
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equation, Eq. (4). We write it with the substitution of Bq by Bq^, as is 
obvious for oblique propagation: 

dv dw Boz dh 

In the gyrofiuid equations [18,19] there are corrections from both E x B 
velocity and an effect of gyro-averaging of fields over the Larmor orbit. 
The E X B drift velocity Ve = c[E x B]/S^ enters the gyrofiuid analogs 
of Eqs. (l)-(3) in the combination V • {u'e^' + v^), where e^/ - is the unit 
sector in z' direction. Upon substituting the electric field from Eq. (35) 
we obtain the following. The first term contributes to the velocity u', 
giving rise to V • {u'egi + v^) = V • v = du/dz, i.e. recovering the usual 
velocity divergence term. The second term does not make a contribution 
because (47r/c)J = VxB = Vxb_L and, thus, J_l = 0. The last two terms 
are both gradient terms. For a plane wave the gradient is reduced to 
V = ez{d/dz). So, both these terms contribute to v^; like [e^ x e^'] ~ e^. 
Upon substituting into the divergence term of Eqs. (l)-(3), it vanishes. 
Thus, we conclude that for a plane single coherent wave, the E x B drift 
does not make a significant contribution to the nonlinear Alfven wave 
dynamics. 

Another effect we should consider is the gyro-averaging of electric 
and magnetic fields over the Larmor orbits of particles. This effect enters 
the equations via a differential operator which approaches unity when 
the ion Larmor radius vanishes. Since it acts on fields only and because 
we construct the model which must reduce to MHD in the case of zero 
ion Larmor radius, finite gyro-radius corrections appear only in Eq. (2), 
leaving Eqs. (1) and (3) unchanged. In the case of oblique propagation, 
Eqs. (l)-(3) are: 

| + ^(P«) = 0, (38) 

du du (dp* d bl + hi Bq^ dh^\ 

Pit. +p^^ = - (Jo) 4r + ^^33^ + - 0, (39) 



dt dz \ dz dz Stt Att dz 

d d\fp*\ Tedq 
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Here we have already used the quasineutrahty condition together with 
Eq. (28) for the electric potential perturbation. This appears on the right- 
hand-side of Eq. (39) as a gradient of p* . This accounts for the charge 
separation in an ion-acoustic wave, and the subsequent gyro-averaging 
of the electric field associated with the electron responce. Eq. (40) is a 
trivial rewriting of Eq. (3) for the quantity p* . 

The operator Jq, which carries out the gyro-averaging operation, is 
a linear operator. It is simply a Bessel function represented in Fourier 
space: 

where V'j^ is the gradient in the plane perpendicular to z' . Going to 
the frame of the wave, it contains a longitudinal projection equal to 
smQ{d/dz) . All other components vanish for a plane wave in a homo- 
geneous plasma. The operator (Jq) is the operator Jq averaged over the 
(Maxwellian) particle distribution function. There are different ways [19] 
to approximate (Jo)- We choose one of the more simple approximations: 

(Jo) - r;/^ c 1 - |v'i + ... - 1 - f sin^ e^, (42) 

where pi — vt/^i is the ion Larmor radius. We neglect, as usual, terms of 
higher order in p?. 

Equations (36)- (40) replace Eqs. (l)-(5) for the case of an obliquely 
propagating wave subject to Landau damping and including finite Larmor 
radius effects. The procedure for derivation of the DNLS-like equation 

describing nonlinear, obliquely propagating Alfven waves is similar to 
that explained in Section II for parallel-propagating waves, but requires 
a bit more algebra. From Eqs. (36) and (37) we finally arrive at the 
DNLS-like equation, Eq. (16), with a new field: 

iP = — ^ — (43) 
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instead of = {bij; + ibiy)/BQ. One can mention that in terms of the field 
the DNLS for obhque waves is more complicated and contains Korteweg- 
de Vries nonlinearities. These terms are proportional to Bq^. and, thus, 
disappear in the case of strictly parallel propagation. 

The velocity perturbation U2 is found from the second-order expan- 
sion of Eqs. (38)- (40). After one integration these equations can be writ- 
ten in the frame moving with the wave as: 



VAP0U2 = P2 + To + ^0 , 

{P2- 



^—ui 1 +X||>C 

VA / " 



V2 «2 



(44) 
(45) 



1 /2 

Upon acting on the second equation with the operator Fq' and substi- 



1/2 

tuting Fq' P2 from the first, we obtain the equation for as follows: 



{(i-/!T;'Vxii£(i-§r;/^ 



7 



^2] 



2 



{i + xii^}rr 



1*1 



(46) 



where /?* = {T^/Ti)j3 and j^^l^ = IV'P — (IV'P)- We solve this equation using 
the inverse operator (22) and the commutation relations (23). Using 

1/2 

approximation for the operator Fq' given by Eq. (42), we finally obtain 
the following equation for U2- 



U2 = VA 



Mi-ATi— +£ M2-iV2 



1*1 



(47) 



where 



Xii 



A2 = 



N2 



(l-/?*) + ^(l-r/7) 
1 ^ (l-2/?-)+xg(l-2/?V7) 

2'^\\ (l-/3*)2+£2(i_^./^)2' 



.7-1 



7 {i-p*)^ + xl{i-P*hr' 



(48a) 

(486) 
(48c) 



19 



Here we used the notation ?7| = (pf/2)sin^0, and Mi and M2 are given by 
Eqs. (34). The solution of this equation can be easily obtained in terms 
of the integral dissipation operator J given by Eq. (17). We thus have: 

«2 = VA |(Mi + M2£) ^ (Ja - J-a) + (iVi + iV2>C) ^1 [1*1'] . (49) 

This is the inhomogeneos solution valid for all values of A^. There is the 
finite homogenious solution when A^ is negative, i.e. /?*<!: 

«2/iom = -"2^ sin(2/|A|), (49') 

where u^in is the constant that can be obtained from initial conditions. 
This solution represents, in a sense, a free modulation wave, i.e. a travel- 
ling modulation of the envelope of the nonlinear Alfven wave train with 
amplitude set by initial conditions. When A^ becomes positive, this so- 
lution of the homogenious equation diverges at infinity like sinh(2;/A). 

The dissipation operator combination can be written as: 



^ y sinh[(^' - z)lh\F{z')dz' , if A2 > ; 
sin {{z' - z)/| A|] F{z')dz', if A2 < . 



-1 



(50) 



It is interesting to note that the kernel of this integral is an antisymmetric 
function of z' — z, as in the resonant particle integral, but is not singular 
at z' = z. Thus, the main contribution to this integral is not from the 
instantaneous position z, unlike the resonant particle operator. It is also 
interesting that this operator does not correspond to pure dissipation, 
because the kernel of the operator Jl itself consists of both symmetric 
and antisymmetric parts. 

When A^ ~ we can write the approximate solution of Eq. (47). 
Using an iterative method, we obtain: 

U2 = VA I (Ml + M2C) + [(A^Mi - iVi) + (A^Ma - N2) ^| • (51) 
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Fig. 3 represents the dependence of A^/?7jj vs. l/j3 for the same (Te/Ti) 
values as in Figs. 2a and 2b. This quantity changes sign at /? ~ /Jo- 
It defines the particular structure of the finite Larmor radius integral 
operator given by Eq. (50). The analogous plots of Ni/rj^^ and A^2/?7| are 
shown in Figs. 4a and 4b respectively. These coefficients are negative 
(there is a region where A^i is positive, but very small), peaked near /?o 
and look similar to graphs of the coefficient M2. 



V. Modulation Instability 



It has been demonstrated [9,26,27] that modulation instability leads 
to the destruction of wave trains and the production of solitons, trans- 
fer of energy from large-scale wave-modes to (damped) small-scale modes, 
the appearance of shock waves and shocklets in the vicinity of bow-shocks 
of planets, etc.. As was already mentioned, waves described by the DNLS 
can be modulationly unstable depending on the sense of their polariza- 
tion. Left-hand polarized waves are unstable when the cubic nonlinearity 
coefficient is positive, and right-hand polarized waves are unstable when 
it is negative. In this Section we investigate the modulation instability of 
dissipative wave packets described by the modified DNLS, Eqs. (16) and 
(12'), and compare the result with that of the DNLS with the resonant 
particle term. We write the general equation as follows: 



d(j) ^ d 
dr dz 



- + Q2JL - (|0P>] )] - = 0, (52) 



where 



and = v\/2Q.i. We consider only left-hand polarized waves (i.e. with 
a minus sign on the third term) which are modulationally stable when 
Qi > in the dissipationless limit, Q2 = 0. 
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Our approach follows that of Mio et.al. [28] and Spangler [14]. Ex- 
pressing the solution of Eq. (52) in the form: 

= ue'"^, (54) 

where ti(r, z) and i?(r, z) are real functions and substituting this into Eq. (52), 
we obtain two equations for the real and imaginary parts, respectively: 

+ Q2U [Jl K - >] ) . + '^f^^zUz + fi^zzU = 0, (55a) 
^ + Qi^zu^ - Qi^z {u^} u + Qi^zuJl {u^ - («^>] - [lUzz + [i^lu = 0, (556) 

where ii = du/dt, Uz = du/dz, and similarly for d and J^. Using the solu- 
tions Uq and ©Oj we superimprose small-amplitude and phase modulation 
such as: 

w = t/o + ew, -i? = Oo + ed, (56) 

where e <^1,Uq = const, Qq = —coq, (©0)2 = ^O; where, in turn, uq and ko are 
the frequency and wavenumber of high-amplitude Alfven wave. Eqs. (55) 
with u and of the form (56) yields (in zeroth order of e) the dispersion 
relation: 

ojQ = iikl- (57) 

Note that no nonlinearity enters the dispersion relation when obtained in 
the frame moving with a wave. 

In first order the linearized Eqs. (55) become: 

t[ + Uz{2QiU^ + 2iiko) + 2Q2U^iJLm)z + ^zzI^Uq = 0, (58a) 
d + 2^zl^ko + 2uQikoUo -Uzzj^ + 2Q2/coC^o = 0, (586) 

To go further, we need to know the Fourier representation of the operator 
Jl- We write: 

e^l'^u{z + QdC, = / e^/^ / dk Uke"'^^+'^^ dC 
-00 J —00 \J —00 / 



00 



dke"''-^^—. (59) 
ik+l/L ^ ' 
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Taking the Fourier transformation of Eqs. (58), with ~ exp{ikz — ilut) , 
we obtain: 

{-ioj + Ar + iAi)uk,u, - Cr'dk,uj = 0, (60a) 

(Br - iBi)uk,co + {~ioJ + iCi)^k,co = 0. (606) 

The coefficients are: 

Ar = 2Q2U^L^k^, Ai = k (2QiU^ + 2fiko + , 

Br = 2Q,Uoko+^-^^ + k -, B, = ^-^^^, 
Cr = k^/iUo, Ci = 2jikkQ. 

Instabihty appears when the imaginary part of the compex frequency u 
is positive. When Q2 vanishes there is no instabihty, so typically 7 <^ uj^- 
Substituting the frequency u = ujr + ^7, we may write an expression for 
the growth rate 7 as: 

IQ2I k^L^ .... 

It is noticable that the growth rate in this case does not depend on the 
amplitude of the nonlinear Alfven wave, unlike the resonant particle case. 
In astrophysical plasma /? ~ 1, so that /c^L^ ^ 1, i.e. the modulation length 
is much smaller the characteristic length of dissipation L. In this case: 

In the opposite case fc^L^ <^ 1: 

VAX\\ko (I-/3/7)' f 

where Im = '^/k is the modulation length. One can see that those waves 
which are modulationally stable in the dissipationless limit become unsta- 
ble when colUsional dissipation increases. The growth rate of modulation 
instability is proportional to the ratio of the coefficients of the nonlocal 
and the cubic nonlinearity terms. When the nonlocal term coefficient is 
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zero, we obtain from the system (60) the known criterion of modulation 
stabihty which depends on the sense of polarization of the wave [10,14]. 

Modulation stability analysis of the modified DNLS with the resonant 
particle integral operator has been implemented in Refs. [14,29]. The 
growth rate for this case is: 



Using Eqs. (34) for the case of low beta, /3 -C 1, we obtain the scaling: 



Modulation instability now appears for waves which are stable in the 
absence of resonant particle effects. The growth rate again is proportional 
to the nonlocal-to-nonlinear term coefficient ratio. Roughly, this growth 
rate is M2/M1 times smaller than the modulation instability growth rate 
predicted by the simple DNLS. 

VI. Conclusions 

In this paper we have considered the influence of kinetic effects on 
nonlinear Alfven wave dynamics. The effects under consideration are 

i) the dissipative (collisional) longitudinal viscous and thermal fluxes, 

ii) resonant particle effects, i.e. Landau damping, iii) the E x B drift 
and finite Larmor radius effects such as those associated with the gyro- 
averaging of fields over the Larmor orbits of particles. All calculations 
have been based on fluid moment and gyrofluid models with Landau 
damping modelled by additional dissipation-like terms. The results ob- 
tained are given below. 

a. Dissipative (collisional) longitudinal viscosity, //y , and thermoconduc- 
tivity, x\\i give rise to an additional, new integral nonlinear term in 




(63) 
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the evolution equation for the nonlinear Alfven waves, and also re- 
solve the (1 — singularity of the derivative nonlinear term in the 
MHD model. The dissipative integral operator is different from the 
integral operator representing resonant particle effects, and is given 
by Eq. (17). 

b. Modulationally stable dissipationless Alfven waves become unstable 
when dissipation is included. The growth rate is proportional to the 
ratio of the coefficients of the integral and nonlinear terms. This 
is similar to the case of the modified DNLS with Landau damping, 
studied in Ref. [14]. 

c. Using the three-moment fluid equations [17], we derived the modified 
DNLS including Landau damping effects. Since in this model dissi- 
pation is no longer an algebraic constant, but an integral operator 
given by Eq. (21), we obtained exactly the same functional form of 
the equation as that obtained from lengthy, full kinetic calculations 
[6,14]. However, the coefficient of both the nonlinear and resonant 
particle operator (nonlocal) terms obtained here are much simpler, 
thus facilitating clear physical interpretation and further analysis. 
Particularly, the phenomenon under study is just the broadening of 
the resonance between acoustic and Alfvenic branches due to the res- 
onant particle interaction with a sound wave. 

d. In the collisionless regime, there is coupling of Alfven and ion-acoustic 
waves, instead of coupling of Alfven and sound waves. The resonant 
particle effect is then strong Landau damping of the ion-acoustic wave 
when Tg/Tj ~ 1. We obtained the temperature ratio dependent coef- 
ficients Ml and M2 of the nonlinear and integral terms, respecitvely. 
These coefficients are much simpler than those obtained full Vlasov 
model [14]. However, they display similar qualitative dependence on 
the /? and Tg/Tj parameters. Figs. 2a, 2b. 

e. In the case of obliquely propagating waves, other kinetic effects may 
be relevant. It was shown that even though the E x B drift veloc- 
ity is not zero (electric field is generated from charge separation in 
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an ion-acoustic wave), it does not contribute to the dynamic wave 
equation. Another effect which does contribute is a gyro-averaging 
of electric and magnetic fields acting on a particle over a Larmor 
orbit of this particle. This effect depends on the angle between the 
wave propagation direction and the ambient magnetic field direction 
as sin^ 0, and vanishes for a strictly parallel-propagating wave. We 
studied this effect using a gyrofluid model constructed from gyrofiuid 
equations [18,19]. This effect results in additional terms in the am- 
plitude evolution equation. These terms are expressed in terms of 
the dissipation integral operator, given in Eq. (17). Hence, this ef- 
fect can be viewed as some specific nonlinear dissipation (similar to 
collisional dissipation) of an obliquely propagating wave due to finite 
Larmor radius. 

We have shown in this paper that MHD models with linear kinetic 
corrections provide correct quantitative description of effects such as the 
dynamics of high-amplitude Alfven waves. It is significant that the ex- 
pressions obtained are, nevertheless, much simpler than those obtained 
from kinetic (Vlasov) calculations, thus facilitating further analysis and 
numerical calculations, and allowing very clear physical interpretation. 
There are some points unresolved in this work. First, when the ampli- 
tude of a wave is high, there are particles trapped in the wave. These 
particles traverse the regions of lower field and bounce between the re- 
gions of higher field. This trapped motion in a beat wave results in am- 
plitude oscillation of this wave at the particle bounce frequency, similar 
to the nonlinear Landau damping process. Obviously, such phenomena 
cannot be represented by a theory based on an assumption that the dis- 
tribution function deviates slightly from a local Maxwellian. Obviously, 
the physical process of trapping will lead to plateau formation and other 
modification of Maxwellian structure. Indeed, waves with b±/BQ ~ 1 will 
have large trapping width and undoubtedly lead to significant distortion 
of the distribution function. It is interesting to note that such distortions 
could potentially mitigate the effects of Landau damping, via local flat- 
tening of (/) at the resonant velocity (for (3 ~ 1). Second, the gyrofluid 
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model of this paper is the simplest possible. It is important to consider 
full set of nonlinear gyrofiuid equations to include gyro-viscosity, finite 
Larmor radius corrections to the Reynolds stress, etc., which are relevant 
for the case of oblique and near-perpendicular propagation. 

As noted above, the kinetic effects which govern the dynamics of a 
collisionless plasma may significantly alter the pictire of nonlinear Alfven 
wave dynamics built upon the MHD plasma model. Inclusion of wave- 
particle resonance effects can significantly alter predictions for modula- 
tional instability (i.e. left vs. right circular polarization dependence of 
growth rate upon parameters such as Te/Tj, etc.). Also, the nonlocal 
structure of the envelope equation which arises from the effects of parallel 
streaming, will likely result in departure from the traditional paradigm of 
collisionless shocks as solitons formed by the competition between non- 
linear steepening and dispersion. In particular, a new time scale, namely 
the ion transit time through the envelope modulation, enters along with 
the steepening and dispersion rates. Strong ion heating will occur, as 
well. Thus, collisionless shock structure may be smoothed or exhibit sec- 
ondary temporal oscillations. These speculations may be easily addressed 
by studies of numerical solution of the (tractable) envelope equations de- 
rived in this paper. The result of these numerical studies will be published 
in Part II of this series. 

Finally, it should be mentioned that an improved understanding of 
nonlinear Alfven dynamics in a compressible plasma may have applica- 
tion in contexts other than collisionless shocks in the solar wind plasma. 
First, a significant fraction of the interstellar medium [30] is hot, col- 
lisionless compressible plasma. Thus, the problems of the galactic dy- 
namo and interstellar turbulence should be approached in the context 
of a collisionless, compressible plasma model. In particular, collisionless 
dissipation via ion heating is a natural mechanism for controlling the 
growth of small scale magnetic energy which has been shown to inhibit 
the mean field dynamo in purely incompressible MHD theories [31,32]. 
Similarly, the processes of wave steepening and shock or soliton forma- 
tion can strongly affect the parallel dynamics of Alfvenic turbulence in 
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the interstellar medium. Such turbulence is thought to be related to 
interstellar scintillations. These issues will be addressed in future publi- 
cations. 
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Figs. 2a, 2b. The same coefRcients as on Fig. 1 in tlie two-fluid model 
vs. 1/P for four T^/Ti temperature ratios (Tg/Ti = 1,3,5,10). 
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Fig. 3. 



Fig. 3. Same dependence as on Figs. 2a, 2b for the normalized quan- 
tity A7?7j. 
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Fig. 4a. 
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Figs. 4a, 4b. Same as on Figs. 2a, 2b for tlie normalized coefficients 
Ni/ril and iVa/ry^. 
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